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稠 油 热 采 三 重 介质 三 渗 模 型 及 压力 动态 分 析 


陈 帅 ,和 孙 可 明 ” 
(1. 辽宁 工程 技术 大 学 力学 与 工程 学 院 ,123000 阜新 ;2. 青 岛 理 工大 学 理学 院 ,266520 青岛 ) 


摘 要 :为 了 研究 碳酸 盐 岩 储 层 稠 油 热 采 的 并 底 压 力 变 化 特征 ,考虑 温度 对 油 度 的 影响 和 地 层 的 应 
力 敏 感性 ,建立 了 三 重 介质 三 渗 数 学 模型 。 数 学 模型 为 非 线 性 变 系数 偏 微分 方程 组 ,采用 全 隐 式 有 
限 差 分 格式 本 用 Newton-Raphson ds 性 差分 方程 组 。 差 分 方程 组 的 Jacobi 短 阵 为 
三 对 角 短 阵 。 利 用 三 对 角 给 阵 的 性 质 求 北 阵 提高 了 和 迭代 速度 。 根 据 计 算 结果 绘制 了 试 并 样板 曲 
线 ， 隆 分 析 。 研 究 结 果 表 明 : ee。 分 为 6 个 阶段 ， 不 同 储 层 中 稠 油 的 夭 度 依赖 
“于 温度 变化 的 敏感 性 不 同 。 具 体 表 现 为 黏度 对 温度 变化 的 敏感 性 越 高 则 储 层 中 压力 的 数值 越 大 ， 
各 窜 流 阶段 出 现 得 越 早 ,而 且 会 缩短 溶洞 向 裂缝 窜 流 进行 的 时 间 并 减弱 其 效果 ;造成 压力 导数 提前 


衰减 。 
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Triple permeability model in triple media and pressure dynamic 
analysis of thermal recovery of heavy oil 


CHEN Shuai' ,SUN Keming” 


(1. School of Mechanics and Engineering ,Liaoning Technical University ,123000 Fuxin, China; 
2. School of Science, Qingdao University of Technology ,266520 Qingdao ,China ) 


Abstract :To study the characteristics of the wellbore pressure changes of thermal recovery of heavy oil in 
carbonate reservoirs ,considering the effect of temperature on viscosity of oil and stress Sensitivity of reser- 
voirs ,the triple permeability mathematical model in triple media has been constructed. The mathematical 
model is of nonlinear partial differential equations with variable coefficients. The difference equations have 
been discretized by the fully implicit finite difference scheme. Newton-Raphson method has been used for 
solving nonlinear difference equations. The Jacobi matrices of the difference equations are tridiagonal ma- 
trix. Using the properties of tridiagonal matrix ,the iteration speed is improved. According to the calculation 
results ,the typical well test curves have been plotted. Sensitivity of parameters has been analyzed. There 
are eight stages of pressure response. The sensitivities of viscosity depend on temperature are different in 
different reservoirs. The higher the sensitivities of viscosity depend on temperature ,the greater the values of 


pressure in reservoirs. When the cross flow stages are earlier ,the time of the cross flow between the frac- 
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tures and vugs is shortened and the strength is weakened,causing the earlier decrease in the pressure de- 


rivative. 


Key words :triple media; carbonate; heavy oil;equation with variable coefficient; well testing 


一 般 碳 酸 盐 岩 储 层 有 裂 颖 和 溶洞 发 育 ,其 渗流 
规律 不 同 于 只 含有 一 种 孔隙 度 的 油气 藏 。 碳 酸 盐 岩 
渗流 流动 规律 比 单一 孔隙 储 层 复杂 ,因为 基质 系统 
(以 下 简称 基质 ) ,裂缝 系统 (以 下 简称 裂缝 ) 和 溶洞 
系统 (以 下 简称 溶洞 ) 间 存在 流体 交换 。 其 流体 流 
动 规律 比 单一 孔隙 储 层 复杂 。 为 此 ,为 真实 反映 碳 
酸 盐 油气 藏 的 储 层 特征 , 需 将 碳酸 盐 岩 简化 为 多 重 
介质 进行 研究 。 文 献 [ 1-2 ] 提出 了 双重 介质 模型 。 
文献 [3] 对 三 重 介质 模型 的 建立 做 了 基础 性 工作 。 
文献 [4] 提 出 了 三 重 孔 阶 介质 概念 。 文 献 [5 ] 指 出 
三 重 介 质 曲线 拟 合 必 须 考虑 溶洞 的 影响 。 文 
献 了 6-9] 对 三 渗 模 型 进行 了 研究 ,并 考虑 了 地 层 应 
} 猎 感性 和 二 次 压力 梯度 的 影响 。 文 献 [10-11 ] 建 
wx 省 求解 了 考虑 启动 压力 梯度 的 分 形 三 重 介质 油 藏 
闫 学 模型 。 文 献 [12] 研 究 了 三 重 介质 地 层 中 油水 
两 相 非 稳 态 窜 流 的 水 平井 模型 。 文 献 [ 13 ] 考虑 压 
化 利 地 层 渗透 率 各 向 异性 ,建立 并 求解 了 斜 井 数学 
模 淹 。 文 献 [14] 建 立 并 求解 了 井 简 与 大 溶洞 连通 ， 
淮 @ 考 虑 了 流动 和 波动 宁 合 的 模型 。 

二 " 随 着 石油 工业 的 发 展 ,发 现 了 许多 碳酸 盐 岩 类 
笛 油 油 藏 , 稠 油 热 采 被 广泛 应 用 。 文 献 [15] 首 先 对 
而 砍 盐 岩 稠 油 热 采 的 多 重 介质 模型 进行 研究 ,把 储 
层 谷 为 蒸汽 波及 区 和 冷 油 区 ,建立 并 求解 了 复合 油 
藏 孙 平 井 试 井 解释 模型 。 但 是 注 蒸汽 改造 后 ,温度 
对 信 体 夭 度 的 影响 从 井 简 到 油 藏 边界 是 连续 变化 
的 ,而 且 基 质 、 裂 多 和 溶洞 的 渗透 性 需要 被 考虑 。 

目前 ,考虑 到 以 上 因素 的 多 重 介质 试 井 模型 还 
没有 被 提出 。 本 研究 考虑 温度 对 流体 秋 度 的 影响 ， 
并 且 考 虑 地 层 应 力 敏感 性 和 二 次 梯度 项 的 影响 , 建 
立 三 孔 三 渗 模 型 ,用 数值 方法 求解 并 对 影响 井 底 压 
力 的 因素 进行 分 析 。 


一 


1 数学 模型 


L 


碳酸 盐 岩 地 层 中 裂缝 和 溶洞 发 育 ,流体 通过 基 
质 、 裂 缝 和 溶洞 流入 井 简 。3 种 流体 系统 间 存 在 流 
体 交 换 。 在 稠 油 热 采 过 程 中 , 莹 汽 的 注入 使 近 井 地 
带 和 和 远 井 地 带 的 温度 产生 差异 。 温 度 的 改变 使 流体 
流动 性 发 生变 化 。 


假设 温度 的 变化 仅 对 笛 油 的 黏度 产生 影响 , 考 
虚 单 层 圆 形 油 藏 中 心 一 口 井 定量 生产 , 因 内 外 边界 
温度 差异 而 引起 储 层 内 的 热传导 。 忽 上 略 因 流体 流动 
产生 的 热 对 流 ; 每 种 介质 的 孔隙 度 随 压力 的 变化 是 
相对 独立 的 ;基质 、 裂 缝 和 溶洞 被 流体 所 饱和 ;流体 
在 3 种 介质 中 的 流动 都 满足 Darcy 定律 。 开 井 前 地 
层 中 各 处 压力 等 于 原始 地 层 压 力 p ;裂缝 和 溶洞 发 
育 良好 ,裂缝 的 渗透 率 大 于 溶洞 的 渗透 率 ;基质 中 的 
流体 向 裂缝 和 溶洞 发 生 帘 流 ,溶洞 中 的 流体 向 裂缝 
发 生 帘 流 , 审 流 是 拟 稳 态 的 。 基 于 这 些 假设 ,本 研究 
考虑 压力 敏感 效应 和 二 次 梯度 项 。 数 学 推导 和 公式 
中 的 物理 量 符号 中 ,下 标 1 代表 裂 颖 ,下 标 2 代表 游 
洞 , 下 标 3 代表 基质 。 三 重 介质 三 渗 模 型 简 图 如 图 1 
所 示 。 


wellbore 


图 1 三 重 介 质 三 渗 模 型 简 图 

Fig.1 Triple permeability flow scheme in triple-porosity media 

假设 注入 蒸汽 时 温度 不 变 , 外 边界 考虑 的 是 储 
层 温度 ,看 作 温度 是 恒定 的 ,内 边界 假设 注 热 井 的 温 
度 是 恒定 的 。 设 井 简 半径 为 7,, 储 层 半径 为 R; 井 简 
处 的 温度 为 7, , 储 层 外 边界 的 温度 为 Th。 根据 稳 态 
热传导 方程 , 储 层 内 的 温度 分 布 满足 下 列 定 解 问题 
及 其 边界 条 件 。 


dd 
r=r, ,7 = 了， ‘by 
Fh 
其 中 :r 为 半径 ,单位 为 m;7 为 温度 ,单位 为 C 。 
解 式 (1) ,得 


S 国 -。 


J 2 
“TIn(R/r.) n(r ri ) ( ) 
定义 理 度 变化 因子 为 
-ld 
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其 中 : 为 流体 黏度 ,单位 为 Pa. s;B 的 单位 为 CC 。 
黏度 变化 因子 表示 流体 黏度 依赖 于 温度 变化 的 
敏感 性 。 
由 式 (3) 得 

MW=MoexpLB(T- Tr)] (4) 
其 中 jw 为 初始 温度 时 的 黏度 ,单位 为 Pa s。 取 BpB 
的 值 为 -0.031 5845% 
将 式 (2) 代 入 式 (4) ,得 


从 -woexp[ 有 | 了 -7T, + mn))] (5) 


考虑 压力 敏感 效应 ,渗透 率 模 数 定义 为 
1 dK,. 
a (6) 
其 中 (i=1,2,3) 为 渗透 率 ,单位 为 m?。 假 设 为 
常数 ,单位 为 Pa 。 则 
人 有 大; = Kexp[ yi(p; 一 po)] ,i=1,2,3 (7) 


其 中 Kio(i=1,2,3) 为 初始 渗透 率 。p,(i=1,2,3) 为 
3 种 介质 中 的 流体 压力 。 


定 流 方程 为 
qd12 = Po(p, -Pp1) 
qi3 = (ps -pi1) (8) 
G23 = op 一 P2) 
内 


其 中 :a 为 裂缝 与 溶洞 之 间 的 量 纲 为 一 的 窜 流 系 
数 ;ai; 为 基质 与 裂缝 之 间 的 量 纲 为 一 的 窜 流 系数 ; 
Qa 为 溶洞 与 基质 之 间 的 量 纲 为 一 的 帘 流 系数 ;gq1,、 
qi3、4w 分 别 为 对 应 的 帘 流 量 ,单位 为 kg/ (ms)。 

将 运动 方程 帘 流 方程 和 状态 方程 代入 连续 性 
方程 (具体 推导 过 程 从 略 ) ,得 到 描述 考虑 温度 影响 
的 三 重 介质 三 活活 流 的 偏 微分 方程 组 , 即 


2 
(po -Pp1) (po -Pp1) 
) | oN (ps -pi) tae "(ps -pi1) 


ot 


2 
(po -p2) (po -p2) 
) | m0) (pp;) + oe m0) (ps —p,) 


(9) 


ot 


2 
(po-p3) (po —p3) 
) | +oaaec 0 Y (pi-p3) +aasec 0 ® (ps -ps) 


2 
OP 1 op op1 
10 2 + + 1 
or 7 Or or 
a =Mopic eel Tv 一 TR + la( /ry) Jryi(ro -Pp1) pi 
$= = U Ww 
5 -2 
9p,» 1 op; Op, 
Kk, 2 + 十 Ya2 | 一 
or r or or 
pacae [mmr et) ps 
Ww 
2 和 
9p3 1 9op; Op3 
30 2 a tyYs 
or r Or or 
CC =Jpsc eel Ty Tat in/rw) ]+73(70 -13) Op3 
和 二 Mo93Ca w 


其 鲜 :p.(i=1,2,3) 为 孔 阶 度 ;c(i =1,2,3) 为 综合 
压缩 系数 ,单位 为 Pa ;i 为 时 间 , 单 位 为 s。 该 模型 
综合 考虑 了 基质 、 裂 颖 和 溶洞 的 渗透 性 ,并 且 考 虑 了 
地 层 的 应 力 敏 感性 和 二 次 梯度 项 的 影响 。 考 虑 温度 
对 流体 黏度 的 影响 ,这 种 影响 从 井 底 到 外 边界 连续 
变化 ,没有 把 储 层 分 区 。 
初始 条 件 为 
Pi(ri| =p(r rsR;i=1,2,3) (10) 
外 边界 条 件 为 定 压 边界 条 件 , 即 
Pi(7,t) |,sR=po(t >0;i=1,2,3) (11) 
考虑 井 简 储 集 和 表皮 因子 的 内 边界 条 件 为 
277,h 0D， 


or 


-0) OP 和 
(Ke po) + Ke™” p0) 十 
r 


3 dp, 
Ken-m 人] | _C 站 =g,(t>0) (12) 


ot 
-» SS_ WM _ 
pw =D1 SIR (13) 


-» SS_W 
pw =P2 STEh (14) 


三 站 = 
Pw =P3 eR (15) 


其 中 :p, 为 井 底 压力 ,单位 为 Pa;h 为 储 层 厚 度 , 单 位 
为 m;g 为 流量 ,单位 为 m/s;C 为 井 简 储 集 系 数 , 单 
位 为 m . s /kg;5 为 表皮 系数 。 

按 表 1 定义 无 因 次 变量 , 式 (9) 化 为 


Opn 1 ppi oppi 
一 2 + Yn | + 
orp 7D orp Orp 


| 


A1e (pp -pop) +Ase’ (pp -Poi) 


0 
Bp [ 90- mir rp |+ ypipp por 


二 CO16 - 
glp ” 
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9 8 pp 
D,| Pm ,1 opp _ 人 | |+ 


Pr 入 
orp rp Orp orp 


YD2PD2 YD2PD2 
Awe (pp -pnm) +A23e (pps -pn ) 


pp[ 0 rp]+ ypapp2 OP D2 
otp ” 
1 pps _ 人 wm) ] 是 
Orp 


= we 


-2 
0 
D, | 下 


orp Tp Orp 


Ase’ (pp —Ppp3) +A, 3 "2 (pp, -pp) 


= ws een[ 6 EC 人 ro]+ynpapns 2 (16) 
tp 


式 (16) 是 非 线性 变 系数 俩 微分 方程 组 ,需要 结 
合 初始 条 件 与 边界 条 件 进 行 数值 求解 。 
表 1 无 因 次 变量 定义 表 


Tab.1 Definitions of dimensionless variables 
符号 
六 膨 扩 表达 式 ( =1,2,3 
i tl , 
到) 2 hio + kao + kao )h 
压 大 PDi 一 一 (po -pi) 
[ Kod 
: 2 hio + hao + hao 
井 底 压力 pup 一 (Po 一) 
Wy Mo9 
A\ (hio + hao + hao )t 
时 间 J - 
> " Mrs Picu + paco + ace) 
径 向 距离 ro 


T Ww 


> My: 

洲 雪 权 类 [7] | PY 

人 9 27h( kio + ko + kao ) 
Ca 


klo + kao + ka0 


也 到 流 系数 A 


2 
2 Q03 7 
孔 缝 窜 流 系数 Al13 一 一 2 
klo + kao + ka0 
Qn37 
洞 缝 帘 流 系 妆 A 一 一 一 一 
洞 缝 窜 流 系数 23 A 
pe Dicti 
储 容 比 Ci ! 


PICu 十 pco 十 p3CcB 


C 
2mhrv (plcu +paco + pacn) 


井 简 储 集 系 数 Co 


kio 
渗 透 率 比 D; 一 -一 一 一 
渗透 率 klo + kao + ks0 
黏度 变化 因子 Bp BTr 
Er T,— Tk 
无 因 次 温度 0 
R 
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初始 条 件 为 


Ppi(rp ,tp) pe =0,(1<r<R/r, ;i=1,2,3) (17) 
外 边界 条 件 为 
Poi(Tpstp) -ar =0,(tp >0;i=1,2,3) (18) 
内 边界 条 件 为 


C dp 一 | D, er pn 4 De -rm PP 平 
D dtp 1 9 2 


rp rp 


De Yops we] 
’ 0 


=1,(t, >0) (19) 


7D rp=1 

0 

Pwp = (mm -rpSe | (20) 
orp rp=1 
0 

Pw = (mm -rpSe | (21) 
QT sd 
0 

Pw = (Pn 一 7DSe “Im | (22) 
orp rp=1 


为 了 消去 二 次 梯度 项 ,对 方程 组 及 其 初始 条 件 
和 边界 条 件 进行 Pedrosa 代 换 -9 , 则 


1 
=- 一 ml[l-yn7(rn 如 23 
本 [1 -ymn(ro,ity)] (23) 


D1 


pm(rp,tp) = 


1 
pmp(rp,tp) = 一 一 Imn[l1-ype(ro)t 如 )] (24) 
Yp2 


pm(rp,tp) = -hl — Yo (rp ,tp)] (25) 


x = Ilnry (26) 
式 中 mE 和 Y 表示 变换 后 的 压力 。 
经 Pedrosa 代 换 后 消去 了 二 次 梯度 项 ,方程 组 变 
为 如 下 形式 。 


2 
ee (1 -yon) -epp[ 0- FT n+ 
e€” Ox” otp 
A Es — Ymé) _ln(1 本 上 
Yp2 YD1 
‘so — yp36) _ln(1 | 
YD1 
20 


pg, (1 yo) pp[ 0 -mir "+ 


otp 


2x 


In(1 2 


_ln(1 | 于 
Yp2 


ln(1 -0 _ln(1 | 
yp2 


otp 


一 Co _ln(1 


| 下 

YD3 

ee 和 | 
Yp3 


| 
al 
peg — yl) el? -Fy e+ 
+ol 
ta 


(27) 
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初始 条 件 为 初始 条 件 离散 为 
nl = 0=6|, .0=0 (28) = = =0 (35) 
外 边界 条 件 为 外 边界 条 件 离散 为 
WN ly=InC Rrs) =é x=InC RIrw) =6 -mt R/r,) =0 (29) uh = = =0 (36) 
内 边界 条 件 为 内 边界 条 件 经 过 文献 [7] 中 的 处 理 离散 为 
dp, Cn ss ) 0 -0 1 
Go (DD tp) ,10) al 
_[_ hl-ymn) on TI” -0 和 -的 -0 
= | yp 5 axl |， -0 | 攻 Ax D, Ax Ax ” 
lIn(1— Cprln(1— g 1 —70 
pw =| - n( 2 - (32) 癌 n( = +] -1=0 (37) 
D1 
1 加 n+l n+l _ gnt+l 
ee =| n( J _ | (33) 2| _ln(1 2 ) 8 2 | _ 
*=0 D2 
i 
2 本 数值 解法 Pt 
(© pf ln(1 ~-ymé) E61 -é0 a 
) 一 阶 导数 采用 向 前 差 商 , 二 阶 导数 采用 中 心 差 | Be |-1=0 38) 
了 对 变 系数 方程 组 进行 全 隐 式 差分 离散 '"" ,得 Ga In ym) SE 全] - 
BD 2 + op) 人 
GS A (Ax)’ n_ on n_ en n+l 本 n+l 
wi | PE oh D, 人 -也 ， 和 -DD, [| 再 
局 一 一 ec nC RIrw I = 
(1 -yom ) Ai 本 n n_ zn 
Le — ymé,; ) _ln(1 -Yo wa | st] -1 =0 2 
NIA， 11 D3 
Bs We ， ee 个 非 线 性 变 系数 差分 方程 与 其 相 
人 应 的 边界 条 件 式 (37) .起 (38) 和 式 (39) 组 成 3 个 广 
人 (Er a + Er ) 程 组 Newton- Raphson 迭代 法 进行 求解 。 
aE 以 式 (34) 中 第 一 式 与 式 (37) 组 成 的 方程 组 为 
n+l n i 入 器 目 . 1 地 
rT CO， | po[。 Fe Ci- De] -5 一 各 性 例 。 ee 本 
(1 -ymé ) 人 fi(nm m1) = 
A 区 一 Ypl771 ) _ln(1 | (40) 
12 
a YDl 二 本 i , fm ,se ,未 号 ) -0 
卫 )3 DA 一 yp 
这 里 及 为 式 (37) 的 左 端 ,A(i =2 ~ ]) 由 式 (34) 
D3 
D, (好 | 20 + ) 中 第 1 式 移 项 得 到 。 用 向量 记号 记 为 
= (Ax)? =) ER’ (41) 
n+l n 一 eo. T 
二 ws pp[ 0- ep -0 + 下 = (A 2 :) (42) 
(1 -ym ) At 则 方程 组 可 以 写 为 
i F(W"' )=0 (43) 
l 7YD2 Yp3 Jacobi 矩阵 为 
到 — yp; ) _ln(1 — Yo30; a (34) J 
- YXYDl YD3 J wn ee Ne (44) 
其 中 :n 从 0 取 到 7T-1;j 从 1 取 到 J 了 -1;J 代表 步 长 2 
的 最 大 值 。 J 
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其 中 地 层 能 量 低 , 井 底 压力 升 高 ,各 帘 流 阶段 出 现 得 越 
(8,0) 早 ,而 且 会 缩短 溶洞 向 询 颖 富 流 进行 的 时 间 并 减 纶 
A 其 效果 ;压力 导数 提前 衰减 。 这 些 影响 会 随 着 稠 油 
Wh 0 对 温度 的 敏感 性 的 增强 而 越 来 越 明显 ,如 图 3 所 示 。 
人 oe 
1 C wD wD "Ip Dp 
J = OL At 。 人 2 ll 6 
m1 Ax At 1 -yom 10' 人 | 
fF ; a | | | | 
Jii1 =- = isti= 一 2( 二 2， 1 pa | 请 | a 
om om e ( Ax) < # ph | 
(i=2,3,.…,] -1;j=i-1), | | | 本 | | 
a D, A tc 
ji = G0 OD Ar) ， 
; ; 1 10 10 10 10° 10 10 10" 
J epp[ 0- 氏 动 富 ] : We 
i 1 一 1 a 汪 、 
的 图 2 并 底 压力 及 压力 导数 曙 线 
By, ( /ee/ ) +1] -CA +As)/( -Yo Uh ) ， Fig.2 The wellbore pressure and pressure-derivative curves 


(© (i1=2,3,.…: 
( 品 经 过 迭代 后 的 新 值 为 
(多 Ee 

MY -TY Fy))® (45) 
其 中 上 为 迭代 次 数 。 

二 计算 过 程 中 每 一 次 迭代 都 要 计算 Jacobi 矩阵 的 
逆 阜 阵 。 由 于 Jacobi 和 矩阵 为 三 对 角 和 矩阵 ,利用 三 对 
秆 所 阵 的 特点 可 以 使 求 逆 矩 阵 更 简单 。 采 用 文献 


[Cj 中 对 三 对 角 和 矩阵 求 道 的 算法 提高 了 计算 速度 。 
3 二 结果 分 析 
© 
.天 描述 温度 影响 下 的 三 重 介质 三 活 渗 流 的 偏 微分 
方 绑 组 是 变 系数 非 线性 方程 组 ,离散 以 后 的 计算 先 
代金 数 比 不 考虑 温度 影响 的 方程 组 迁 代 次 数 多 ,这 
种 现象 随 着 5 绝对 值 的 增 大 而 越 来 越 明显 。 

在 双 对 数 坐 标 中 画 出 井 底 压 力 及 压力 导数 曲 
线 ,如 图 2 所 示 。 可 以 划分 为 6 个 流动 阶段 :阶段 1 
为 并 筒 储存 阶段 ,压力 和 压力 导数 曲线 呈 单 位 斜率 
直线 ;阶段 2 为 表皮 因子 影响 阶段 ,压力 导数 曲线 呈 
驼峰 状 ;阶段 3 为 溶洞 向 裂缝 窜 流 阶段 ,压力 导数 曲 
线 出 现 止 陷 ;阶段 4 为 基质 向 裂缝 定 流 阶段 ,压力 导 
数 曲 线 出 现 凹陷 ;阶段 5 为 基质 向 溶洞 窜 流 阶段 , 故 
力 导数 曲线 出 现 四 陷 ;阶段 6 为 稳定 流动 阶段 ,压力 
导数 急剧 下 降 , 并 底 压力 图 线 为 水 平 直 线 。 

对 计算 结果 进行 参数 敏感 性 分 析 如 下 。 

1) 由 于 不 同 储 层 中 稠 油 对 温度 变化 的 敏感 性 不 
同 ,B, 绝 对 值 也 有 所 不 同 。 注 燕 汽 对 储 层 进行 改造 
后 , 注 燕 汽 后 流体 的 流动 性 增强 ,定量 生产 时 消耗 的 
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图 3 Bb 对 无 因 次 井 底 压力 及 压力 导数 的 影响 


Fig.3 The effect of B, on dimensionless wellbore pressure 


and pressure-derivative 
2 ) 渗透 率 模 数 越 大 ,渗透 率 变化 对 压力 改变 越 
敏感 。ym 越 大 , 井 底 压 力 越 大 ,压力 导数 的 数值 也 
越 大 ,如 图 4 所 示 。 


10 Pp:— Yo-0.1 — yw=0.15 一 一 y=0.2 
“yoi=0.2 


Case Yo 0. we yor0.15 oo 


图 4 yoi 对 无 因 次 井 底 压力 及 压力 导数 的 影响 


Fig.4 The effect of ym on dimensionless wellbore pressure 


and pressure-derivative 

3) 窜 流 系数 和 ;影响 各 审 流 阶 段 发 生 的 时 间 , 其 
值 越 大 ,溶洞 向 肌 缝 窜 流 阶段 进行 的 时 间 越 短 ,基质 
向 裂缝 究 流 阶段 进行 的 时 间 越 长 而 且 效 果 越 明显 ， 
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如 图 5 所 示 。 6) 表 皮 系 数 $ 越 大 , 井 底 压力 越 大 ;表皮 因子 影 
必 Pio:—hs0.1 一 -015 一 -2 一 0.2 啊 阶 段 持续 的 时 间 越 长 而 且 压 力 导 数值 越 大 ; 游 洞 
| 向 裂缝 罕 流 阶段 出 现 得 越 晚 ,持续 时 间 越 短 而 且 效 
5 果 越 弱 , 如 图 8 所 示 。 
i eh WwW 10? 0 
~ AAA SS Pup” to/Co:  S=1 
3 
1 
本 
es < 


图 5 As 对 无 因 次 井 底 压 力 及 压力 导数 的 影响 
Fig.5 The effect of A on dimensionless wellbore pressure 
and pressure-derivative 

4) 储 容 比 w; 影 响 各 帘 流 阶段 出 现 的 早晚 和 持续 

时 间 , 图 6 给 出 了 w, 对 曲线 形态 的 影响 ,w, 不 变 ,w 越 

小 ; 越 大 ,溶洞 向 裂缝 窜 流 阶段 持续 时 间 越 长 ; 基 
质 则 裂缝 窗 流 阶段 出 现 得 越 史 ,持续 时 间 越 短 。 
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图 6 w; 对 无 因 次 井 底 压 力 及 压力 导数 的 影响 


FE 下 .6 The effect of w, on dimensionless wellbore pressure 


naXiIv 


and pressure-derivative 

- 己 5) 间 简 储 集 系 数 Cs 对 早期 曲线 的 形态 影响 比 
较 只 ,Ct 越 大 ,早期 的 井 底 压 力 越 小 , 井 简 储 存 阶段 
持续 时 间 越 长 ;溶洞 向 裂缝 富 流 阶段 出 现 得 越 晚 , 持 
续 时 间 越 短 而 且 窜 流 效 果 越 弱 ; 基质 向 裂缝 窜 流 阶 
段 出 现 得 越 晚 ,持续 时 间 越 短 , 如 图 7 所 示 。 


Per; 10 
Pw " th/Corem C10 一 


i/ 


图 7 Cb 对 无 因 次 井 底 压力 及 压力 导数 的 影响 


Fig.7 The effect of Ch on dimensionless wellbore pressure 


and pressure-derivative 


R/C 
图 8 5S 对 无 因 次 井 底 压力 及 压力 导数 的 影响 


Fig.8 The effect of S on dimensionless wellbore pressure 


and pressure-derivative 
* 
八 
4 结 论 


1) 碳 酸 盐 岩 油 藏 考虑 为 三 渗 三 重 介质 建立 试 井 
模型 ,压力 响应 分 为 6 个 阶段 : 井 简 储存 阶段 ,表皮 
影响 阶段 ,溶洞 向 裂缝 帘 流 阶段 ,基质 向 裂缝 窜 
流 阶 段 ,基质 向 溶洞 富 流 阶段 和 稳定 流动 阶段 。 

2) 各 个 参数 对 压力 响应 曲线 的 各 阶段 有 不 同 的 
影响 。 注 蒸汽 对 储 层 进 行 改造 后 , 注 燕 汽 后 流体 的 
流动 性 增强 ,定量 生产 时 消耗 的 地 层 能 量 低 ; 井 底 压 
力 升 高 ;各 窜 流 阶段 出 现 得 越 早 ,而 且 会 缩短 溶洞 向 
裂缝 帘 流 进行 的 时 间 并 减弱 其 效果 ;压力 导数 提前 

3) 不 同 储存 中 的 稠 油 恭 度 对 温度 的 敏感 性 不 
同 ,这 种 敏感 性 用 Bn 来 刻画 ,其 绝对 值 越 大 , 则 黏度 
对 温度 的 变化 越 敏感 ,Bb 对 井 底 压 力 及 压力 导数 影 
响 就 越 明 显 , 差分 离散 后 计算 时 迭代 的 次 数 也 


会 增多 。 
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